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ABSTRACT 

Phased array feeds (PAFs) for reflector antennas offer the potential for increased reflector field of view 
and faster survey speeds. To address some of the development challenges that remain for scientifically 
useful PAFs, including calibration and beamforming algorithms, sensitivity optimization, and demonstra- 
tion of wide field of view imaging, we report experimental results from a 19 element room temperature 
L-band PAF mounted on the Green Bank 20-Meter Telescope. Formed beams achieved an aperture ef- 
ficiency of 69% and system noise temperature of 66 K. Radio camera images of several sky regions are 
presented. We investigate the noise performance and sensitivity of the system as a function of elevation 
angle with statistically optimal beamforming and demonstrate cancelation of radio frequency interference 
sources with adaptive spatial filtering. 



1. Introduction 

The generation of radio astronomy instruments now 
coming on line will fill much of the observational 
parameter space — frequency coverage, instantaneous 
bandwidth, system noise temperature, angular resolu- 
tion, and time and frequency resolution. The two re- 
maining open-ended parameters where fundamentally 
new science will be explored are collecting area and 
field of view. More collecting area will allow us to ob- 
serve known phenomena much deeper in the universe, 
and greater fields of view open the possibility of find- 
ing new phenomena, such as transient radio sources 
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and rare types of pulsars, by making large-area sky 
surveys much more efficient. Several major instru- 
ment development projects are underway worldwide 
to develop radio cameras, variously referred to as ac- 
tive, phased, beamforming, or smart arrays to distin- 
guish them from the more conventional independent- 
pixel feed-horn arrays which sample less than l/16th 
of the available sky area within the array's field-of- 
view. These instruments will utilize focal plane phased 
array feeds (PAF) which can electronically synthe- 
size multiple, simultaneous far field beams for com- 
plete coverage of the field of view without loss of 
sensitivity in each beam. PAFs are in development 
or in planning stages for single dish instruments in- 
cluding the Green Bank Telescope (GBT) as well as 
the Westerbork Synthesis Radio Telescope (WSRT) 
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dVerheiien et al J2008h and the Australian Square Kilo- 
meter Arra^jASJCAP}JPathfinder synthesis imaging 
array ds. G. Hay et al.ll2007h . 

As the cost of signal processing for beam-forming 
and the correlation of multiple beams in aperture syn- 
thesis arrays allows, wider fields of view for SKA and 
its predecessors will open up new science, such as sur- 
veys of neutral hydrogen at high redshifts to study the 
evolutionary history of the universe. First science that 
will be enabled by PAFs with modest signal processing 
bandwidths include observations of the trajectories and 
abundance of high- velocity H I clouds interacting with 
the d i sk of the Milky Way dWakker & van Woerdenl 
1997tlPutmanll2006tEIockman et al.ll2008l) . surveys of 
gas clouds outside other galaxies covering large ar- 



eas dD. A. Thilker et al.|[2004t iBraun & Thilkei 



2004 



M. Grossi et alj 2008 ). and studies of the kinematics 



of extended HI clouds contained by galactic groups 



dS. Borthakur et al.ll2008tlBraun & Thilkerl2004l;letal 
2008al) . These types of observations are time consum- 
ing with single-pixel instruments and could be con- 
ducted more rapidly and to greater depth using a PAF 
with wide field of view. 

Demand for use of the GBT to study the inner parts 
of the Milky Way disk is extremely high. Wider fields 
of view would enable additional studies of the galac- 
tic center and objects such as the physical conditions 
and structure of the recen tly discovered "far" coun - 
terpart to the 3-kpc arm (IDame & Thaddeusl 120081) . 
A PAF on the GBT would be of immediate use in 
analyzing the complex chemistry and relationship to 
dust evolution of the diffuse interstellar medium us- 
ing 18cm emission from OH, and in understanding 
the medium to minimize the effect of foreground mat 
ter on cosmological inv e stigati o ns dBlack & Dalgarn 
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et al.ll2008bh . 

As compared to single-pixel instruments, PAF 
beam formation requires a substantial amount of sig- 
nal processing for array calibration, operational beam- 
forming, and image formation. Calibration and beam- 
forming algorithms and per formance criteria f or ar- 
ray feeds were surveyed by Ijeffs et alj d2008bl) . Ra- 
dio frequency interference (RFI) mitigation using 
adaptive nulling algorithms was demonstrated by 
Nagel et alj (12007b . Other recent results demonstrat- 
ing feasibility of PAF-based interfe rence mitigation 
in astronomical applications include djeffs & Warnick 
20081 120071 Ijeffs et al.ll2008bl: lLandon et al.ll2008allhb 



Warni ck et al I l2007h 




Fig. 1. — 19 element single polarized PAF and front 
end box containing analog downconverters and re- 
ceiver front ends, mounted on on Green Bank 20- 
Meter Telescope (October, 2007). 



In this paper, we present results from an L-band 
prototype PAF on the Green Bank 20-Meter Tele- 
scope. Early radio ca mera images obtained with the 
PAF were reported by I Warnick et alj (|2008). A more 
detailed study of performance metrics including sensi- 
tivity and system noise are given here. A major con- 
cern with PAF development has been optimizing the 
interfac e between the array and front end receiver elec - 
tronics dWarnick & Jensenll2007l:llvashina et alll2008l) . 
Due to mutual coupling between array antenna ele- 
ments, standard impedance matching techniques used 
for single-port antennas must be extended to multi- 
port systems using the theory of active impedance s 



dE. E. M. Woestenburgl I200H IWarnick et afll2009bl) . 



The performance improvement that could be realized 
using these methods are explored. 

PAFs offer the potential for optimizing beam pat- 
terns to improve sensitivity as the noise environment 
changes or to maximize system performance for a 
given type of observation. We study the dependence 
of sensitivity and efficiency of formed beams on the 
reflector tipping angle as sky noise increases near the 
horizon. Results of RFI mitigation experiments are 
also shown, in order to demonstrate the potential for 
operational use of RFI nulling algorithms in observa- 
tions. These results are strong evidence for the feasi- 
bility of a high sensitivity, wide field of view PAF. 

2. Experiment Description 

The characterization techniques, calibration meth- 
ods, and beamforming algorithms used in the reported 
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Fig. 2. — Ground shield and PAF in sky noise mea- 
surement facility with absorber being lowered over the 
array for Ti wt measurement measurement (July, 2008). 



I- 



Fig. 3. — Dipole element diagram. 



work were developed with a prototype 19 element 
array (Fig. [TJ of single polarized, thickened dipoles 
backed by a ground plane, with a 20 channel non- 
cryocooled receiver and data acquisition system for 
real time sampling and streaming to disk. The array 
was initially characterized using the BYU "Very Small 
Array" 3-meter prototype platform and the NRAO an- 
tenna range in Green Bank, WV. In October 2007 
and July 2008 the prototype PAF was mounted on the 
Green Bank 20-Meter Telescope to measure aperture 
and spillover efficiency, demonstrate multiple beam 
formation, and test RFI mitigation algorithms. 



2.1. Array Feed and Data Acquisition 

The PAF was located at the focus of the Green Bank 
20-Meter 0.43 f/D reflector. Array elements were 
balun-fed half-wave dipoles tuned to 1600MHz A/4 
above a ground plane. The dipoles achieved a 570 
MHz impedance bandwidth (-10 dB reflection coeffi- 
cient) or a 35% bandwidth relative to the center fre- 
quency, which is modest but adequate for initial tests. 
The elements are arranged in two concentric hexag- 
onal rings around a center dipole. The hexagonal 
grid permits 0.6A spacing between dipoles, which is 
slightly farther apart than the typical A/2 for a rect- 
angular grid. This increased inter-element spacing re- 
duced mutual coupling without undersampling or pro- 
ducing grating lobes. A diagram of the dipole element 
is shown in Figure [3] Additional details were given by 
Nagel etal.ld2007l) . 

The dipoles were connected through the array 
ground plane to uncooled LNAs (Ciao Wireless, Inc., 
Camarillo, CA) with 33 K noise temperature and 41 
dB gain at 1600 MHz. The LNAs were measured by 
L. Belostotski (U. Calgary) to have the following noise 
parameters at 1600 MHz: T min = 33 K, R a = 3.417, 
and r opt = 0.07Z90 . 

Receivers in the front-end box at the dish focus 
performed a two-stage downconversion to a final IF 
centered at 2.8125 MHz with a 3 dB bandwidth of 
425 kHz. IF signals were fed by coaxial cable into 
the pedestal room of the telescope to a data acquisi- 
tion system in a shielded rack. IF voltages were syn- 
chronously baseband-subsampled using 12-bit quanti- 
zation for high dynamic range at a sample rate of 1 .25 
Msamples/sec per channel for 20 channels, using five 
commercial data acquisition cards with four channels 
per card. A single network server class PC housed 
all acquisition cards and a high speed array of hard 
drives for real-time data streaming to disk. Beams 
were formed in post-processing on the raw voltage 
samples. 

2.2. Isotropic Noise Measurement Facility 

Measurements of the array isotropic noise re- 
sponse can be used to determine the aperture effi- 
ciency and sys tem noise temperature associated with 
formed beams dWarnick & Jeffsll2008b . The array was 
mounted in a facility with a retractable roof to pro- 
duce a noise field that is approximately isotropic at 
the temperature of cold sky over the array's significant 
response directions. A copper screen connected to the 
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Interfering 




High gain primary antenna 



Fig. 4. — The primary advantage of FPA telescopes is 
increased field of view provided by multiple, simulta- 
neously formed beams. Spatial cancelation of interfer- 
ing signals is also possible. 

array ground plane at a slope of 30° provided partial 
shielding of the array from noise radiated by ground, 
horizon, and nearby buildings. This provided the cold 
load for a F-factor array isotropic noise temperature 
measurement. For the hot load, a 2.4m x 2.4m sheet 
of RF absorber was lowered over the copper screen 
(Fig. [2). The computed array output voltage correla- 
tion matrices for the two configurations were used to 
characterize the system noise for the results given in 
Section O 

3. PAF Beamforming 

In this section we discuss principles and method- 
ologies for practical beamforming with an astronomi- 
cal phased array feed instrument. Figure [4] illustrates 
how digital beamforming techniques enable a PAF to 
form multiple simultaneous beams to increase the in- 
stantaneous field of view for the telescope. 

3.1. Signal Model 

We assume that the array output signals are pro- 
cessed in narrow subbands such that B D/c, where 
B is bandwidth, D is the PAF diameter, and c is the 
speed of light. After downconversion and sampling, 
the M element PAF produces a length M x 1 complex 



baseband data vector at time sample n: 

x[n] = as[n]+n[n] (1) 

and where s[n] is the signal of interest (SOI) and n[n] 
is the total system noise vector seen at the array. In the 
presence of interference, n[?i] will also include com- 
ponents due to a man-made detrimental signal d[n] as 
illustrated in Fig. 0] Vector a is the normalized ar- 
ray response to a unit amplitude point source in the far 
field direction corresponding to s[n]. 

Assuming zero mean wide-sense stationary signals 
and statistical independence between distinct signal 
and noise sources the array covariance is 

R = E{x[n}x H [n}} = R s + R„ (2) 
R s = crfaa H 

Rn = Rrcc + Rsp + Rsky + Rloss + Rint (3) 

where is SOI signal power (variance), and R lcc , 
R sp , R s ky, Rioss and R; nt are receiver, spillover, sky, 
antenna resistive loss, and interference noise covari- 
ance matrices, respectively. R s ky refers to non-SOI 
noise sources in the beam main lobe, including at- 
mospheric, cosmic microwave background (CMB) and 
galactic background (GB) radiation. Models for each 
of these terms, including their spatial structure across 
the array and effective noise temperatures seen at the 
beamformer output, will be introduced in Section [4] 
Except in Section [7] where we consider interference 
canceling, we will assume R; nt = 0. 

R is estimated for the jth short term integration 
(STI) window with sample covariance matrix 

R« = - £ x[n]x H [n] = -X 3 Xf (4) 

n=jL 

X, = [ x [jL],x[jL + l],--- ,x[(i + l)i-l]]. 

where Xj is the L sample long STI data window. The 
required length L depends on the signal scenario and 
desired operational mode. If adaptive cancelation is 
used to suppress moving interference as discussed in 
Section|7] beamformer weights (computed from RW) 
must be updated rapidly with L short enough that Ri n t 
does not change over the STI window. A slower update 
rate, and thus longer L, is used to adapt to spillover 
and sky noise due to pointing changes, as discussed 
in Section [6] Finally, in stable signal conditions, L 
can be very long (on the order of minutes or hours) to 
reduce sample estimation error and yield accurate cal- 
ibration or high sensitivity beamformer solutions. In 
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these cases superscript (j) will be dropped for sim- 
plicity when this does not lead to an ambiguous inter- 
pretation. 

3.2. Calibration 

Since multiple simultaneous beams are formed with 
a PAF as shown in Fig. [4] a calibration for the sig- 
nal array response vector a^ must be performed for 
each direction, fij, corresponding to each formed 
beam's boresight direction, and any additional direc- 
tions where point constraints in the beam pattern re- 
sponse will be placed. Periodic re-calibration may be 
necessary due to strict beam pattern stability require- 
ments, to correct for differential electronic phase and 
gain drift, and to characterize changes in receiver noise 
temperatures. 



We proposed a calibration procedure in (IJeffs et al 



2008b). which is improved here with a lower error, 
noise whitening approach in step [3] 
Calibration Procedure: 

1. Noise covariance, R n : Steer the dish to a rela- 
tively empty patch of sky so x[n] n[n], and 
collect a long term (large L, e.g. 10 minutes) 
sample covariance estimate for the noise field 
R n using (0). 

2. Signal-plus-noise covariance s,Tli: While track- 
ing the brightest available calibration point 
source, steer the dish to calibration angle fli 
(relative to this source). The observed signal 
model is x[n] = ais[n] + n[n], where is the 
desired calibration vector at direction 51^. Cal- 
culate Ri using dUi and the same L as in step 

m 

3. Array response, a^: Compute a.; = R n u max , 
where u max is the dominant solution to the gen- 
eralized eigenequationRiU max = A max R n u max . 
This noise-whitening method produces a lower 
variance estimate with less bias arising from 
correlated noise in R n (caused by inter-e lement 
mutual coupling) than does the method o flJeffs et al 
d2008bl) . 



4. SOI covariance, R S) j: Compute R s .; = A m£ 

5. Form calibration grid: Repeat steps [2}]4] in a 
grid pattern corresponding to the desired dis- 
tribution of beam centers and constraint points, 
e.g. for spherical angles {fl; 1 < i < K}. 



3.3. Beamformer Design 

In order to form beams with a PAF, a design proce- 
dure is required to obtain a set of beamformer weight 

(i) 

vectors w - , from which the output for the ith beam 
is computed as 



Di[ n \ 



U) H r I 

w; x ?i 



< i < K, j 



n 

LJ 



(5) 



A distinct beamformer weight vector w 4 is used for 
each main lobe steering angle Slj. For beamformer 
weights that are fixed during each STI, the dependence 
on the STI index j can be dropped. In a practical PAF 
observing scenario the beams are steered in a rectan- 
gular or hexagonal grid pattern with crossover points 
at the -1 to -3 dB levels. Proper design of a set of 
beamformer weights allows one to steer the main 
beams, control beam shape and sidelobe levels, opti- 
mize sensitivity, and direct placement of nulls towards 
interferers. A thorough introduction to array signal 
processing and beamformer design can be found in 
ilVan Veen & Bucklevll 19881: IVan Treesll2002l) . 

In general, beamformer weights can be designed us- 
ing an a priori electrical model for the array, character- 
ization data measured before installation of the array, 
or in situ calibration data obtained with a procedure 
such as that outlined in the previous section. As will 
be argued below, only the latter of these beamformer 
design procedures is likely to be viable for astronomi- 
cal PAFs. 

Relative to phased arrays for wireless, communi- 
cations, radar, sonar, and other applications, an im- 
portant distinguishing characteristic of astronomical 
PAF beamforming is that given the current state of the 
art the instrument performance with data independent 
pre-computed beamformer weights is inadequate. For 
simulation studies, the conjugate field match beam- 
former has been commonly used to obtain beamformer 
weights, but the only demonstrated experimental ra- 
dio camera images formed to date use statistically 
optimal, adaptive, data dependent beamforming al- 
gorithms (Verheiien et ai1l2008l: lOosterloo et alj|2008t 



Jeffs et al.l2008a|jwarnick et al.l2008l:ljeffs et al.l2009l) . 
The primary reason for this is that the perturbations in 
array element patterns, complex receiver gains, and ca- 
ble phase lengths that occur when the array is mounted 
on the reflector mean that numerical models or antenna 
range characterization are not accurate enough to al- 
low precomputation of beamformer weights for high 
sensitivity or controlled pattern shapes. 
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In applications with less stringent sensitivity and 
stability requirements, it is possible to design 2-D 
array beamformers deterministically to match some 
fixed array response criterion. A number of classical 
methods are used, including windowed beamforming 
and numerical response o ptimization over a dense grid 
of far field samp le points ( Van Veen & Buckley! 1988 ; 
Van TreeslExm It has been shown in simulation that 



using dense calibration response and numerical opti- 
mization it is possible to design well behaved astro- 
nomical PAF beam shapes with little coma, distortion, 
or polarizatio n rotation as beams are steered off-axis 
(IWilUsll2009h . All such methods require an accurately 
calibrated array (by antenna range measurements or 
mathematical modeling) over the entire response field 
of interest. 

With astronomical PAF instruments it is not possi- 
ble a priori to obtain the full-sphere calibration data 
required to design beamformer weights for a given re- 
sponse criterion. Antenna range measurements of the 
bare PAF element response patterns do not account 
for such effects as interactions with support structures. 
Numerical modeling of the array and dish combina- 
tion provides an excellent qualitative match, but our 
studies of precomputed beamformer weights indicate 
that modeled results differ sufficiently in fine phase 
and gain detail over the field of view. Thus beam- 
former weights designed using numerical electromag- 
netic simulations lead to on-dish real data beamformer 
results with lowered sensitivity and beampatterns that 
deviate significantly from the desired shape. 

Detailed modeling of the array elements, back- 
plane, array structure, and dish reflector system can 
provide excellent qualitative match and prediction 
of achievable real-world performance when used to 
compute beampatterns in simulation. The fine detail 
mismatch though precludes use of these beamformer 
weight designs directly on the real system. We specu- 
late that the discrepancies are in subtle fine scale dif- 
ferences between physical geometry and the model. 
The simulations can predict how an array and dish 
of this design will perform, but cannot match the exact 
phase and gain variation subtleties of a given telescope 
necessary to design a controlled beampattern. Details 
of reflector surface roughness, support structure scat- 
tering, small variations in element position and con- 
struction, noise temperature variation between LNAs, 
etc. can only be captured statistically in the models, 
but the exact unknown values cannot be matched to 
the real values for a specific dish and array. 



Since a priori calibration is insufficiently accurate, 
in situ calibration is required. Even with bright cal- 
ibrator sources, due to limited SNR and integration 
time requirements it is not practical to obtain calibra- 
tion vectors a, with adequate density or coverage be- 
yond the first one or two sidelobes (although correla- 
tion of the array outputs with a second high sensitiv- 
ity antenna could be used to improve coverage). Er- 
rors introduced in transferring a beamformer design 
from simulation invariably lead to unacceptable reduc- 
tion in real-world sensitivity. Further, fine control of 
beam pattern shape to reduce temporal variation of 
the beam response ("pattern rumble") and insure uni- 
formity across all pixels cannot be attained with pre- 
computed data independent designs. We have investi- 
gated these deterministic design methods, but perfor- 
mance using real experimental data is poor. We are 
currently studying methods for translating a desirable 
beamformer design from a detailed numerical simula- 
tion to the corresponding actual PAF array using avail- 
able in situ calibrations, but we have found that simply 
correcting for the per-channel phase and gain calibra- 
tion is insufficient due to the effects of array mutual 
coupling. 

For these reasons, rather than designing the beam 
pattern shape over a full sphere using classical beam 
pattern fitting methods, we have proposed using sta- 
tistically optimal, data-dependent beamfo rming de 



sign m ethods for PAF array beamforming (Je ffs et al 



l2008bh . This approach uses calibration grid data over 
a limited field of view near the reflector boresight to- 
gether with the measured noise response of the array 
to design beamformer weight coefficients. In effect, 
the calibration grid data controls the pattern main lobe 
for each beam, and the measured array noise response 
is used to optimize the deep sidelobes to minimize 
spillover noise. Since receiver noise for a PAF is cor- 
related, the statistically optimal beamformer design 
approach also suppresses receiver noise to a degree as 
well. 

We have used two of the several well known statis- 
tically optimal beamforming algorithms: the linearly 

constrained minimum variance (LCMV) and maxi- 

mum sensitivity (max SNR) beamformers (iVan Veen & Buckley 
ll988llVan Treesll2002l) . LCMV minimizes total output 
power while satisfying a set of linear constraints. The 
LCMV optimization problem is given by 

wlcmv.i = argminw Rw subject to w H Ci = f 

w 

(6) 
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For example, using Gj = [e 



, a.j p ] and f 



[1, n, • ■ • , rp] T constrains the beam main lobe to have 
a response of unity in direction fij and r p in P other 
directions, flj . The single constraint minimum vari- 
ance distortionless response (MVDR) beamformer is a 
special case of LCMV obtained when Cj = a^, and 
f = 1. 

Using calibration data only, © can be solved ap- 
proximately as 



WLCMV,! = R n Cj[C i R n Cj] f 



(7) 



We refer to this solution as a "fixed-adaptive" beam- 
former since though it is statistically optimal for the 
calibration noise environment, it is computed using a 
noise covariance measured before or after the observa- 
tion phase and and does not use the current array co- 
variance during the observation. For noise field track- 
ing or moving interference canceling we use a fully 
adaptive mode, replacing R n in Eq. (0 with RW and 
updating the weight computation for each STL 

The maximum sensitivity beamformer is defined as 



w H R s w 



w m sNR,i = arg max — jt 



w w^R n w 



(8) 



The maximization in Eq. (0 invokes the generalized 
eigenvector problem 



Rs,iW m gNR,i — A max R n W 



mSNR.i 



(9) 



for which the solution using only calibration data 
yields a practical fixed-adaptive beamformer. Rapid 
tracking to adapt to noise field evolution can be ac- 
complished by replacing calibration R n in Eq. ([8]) 
with periodically updated "off-source-steered" esti- 
mates R ft obtained during observation. This method 
is not well suited to interference mitigation since 
off-steering changes the spatial structure of R; nt in 
R n . The maximum sensitivity beamformer has been 
used for most experimental PAF observations reported 
to date jWarnick et al.1 120081: IVerheiien et all l2008t 
Oosterloo et al J l2008l:IJeffs et alJl2008all200^r 

Some of the PAF science applications described in 
Section[T]require beams with very low sidelobes. The 
maximum sensitivity beamformer ^ does not result in 
lowest possible sidelobes. In principle, the same cali- 
bration grids used to generate the signal response cor- 
relation matrices R s $ can also be used to design beams 
with controlled sidelobes. By definition, these beams 
will result in lower sensitivity than (|8), but the over- 
all instrument performance should be better for obser- 



vations requiring low sidelobes. We are currently in- 
vestigating beamformer design procedures that incor- 
porate beam shape metrics while still maintaining as 
high a sensitivity as possible. The flexibility to opti- 
mize beams for different science applications in post- 
processing represents an advantage of PAFs over fixed 
feeds. Sidelobe pattern "rumble" or variation between 
dishes is of particular concern in PAF synthesis imag- 
ing as is planned for the ASKAP array and PAF up- 
grade to the WSRT Stable know sidelobe patterns per- 
mit high dynamic range imaging. Since calibration 
grids are only possible (due to SNR limitations) out 
to the first sidelobe or two, it will not be possible to 
strictly control deep sidelobe patters with the beam- 
former weights. On the other hand, outside the first 
few sidelobe rings, the highly attenuating pattern of the 
dish aperture itself dominates over anything the array 
feed can control. This is a positive promising aspect 
of sidelobe pattern rumble control that needs further 
study. 

4. Performance Metrics and Noise Models 

4.1. Sensitivity and Efficiency 

For a phased array, sensitivity, efficiencies, and sys- 
tem noise temperature are beam-dependent. The pri- 
mary figure of merit for a formed beam in the absence 
of interference is the beam sensitivity 



A 
71, 



cir 



2k h 



10-26_p 



-SNR = 



2fch 



w H R s w 



10-26 F s w H K n w 



(10) 



where A e s is the effective receiving area for the PAF 
illuminated dish using beamformer weight w, T sys is 
the beam equivalent noise temperature, k\, is Boltz- 
mann's constant, and F s (Jy) is the flux density of a 
signal source of interest. Sensitivity is related to radia- 
tion, aperture, and spillover efficiencies (T/ ra d, ??ap, and 
T] sp respectively) according to 



A 



eff 



^rad^ap^p 



Tsys ^?rad(7 n sp ~t~ Psky) ~l~ (1 ??rad)7p T Tec 

(ID 

where A p is the physical aperture area; T sp is the 
equivalent spillover noise temperature; T s k y = T cm b + 
T g b + T atm is the combined non-SOI noise tempera- 
ture in the beam main lobe due to cosmic microwave 
background radiation, galactic background, and atmo- 
spheric noise; T rec is the beam equivalent receiver 
noise temperature; and T p is the physical temperature 
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Fig. 5. — The sources of noise that contribute to the 
signal received by a radio telescope include spillover 
noise originating from outside the edge of the reflector 
dish, atmospheric noise collected by the main lobe of 
the antenna radiation pattern, and receiver noise. The 
equivalent receiver noise temperature includes the ef- 
fects of mutual coupling between antenna array ele- 
ments and impedance mismatches between the array 
and receivers. 



of the array antenna. These constituent noise tempera- 
tures arise from the corresponding array noise covari- 
ance terms R sp , R s k y , Riossi an d Rrcc- 

Defining antenna figures of merit is straightforward 
for a passive antenna, but for an active receiving ar- 
ray the standard definitions for antenna terms are not 
directly applicable. The IEEE standard definition of 
ape rture efficiency can be e xtended to active arrays us- 
ing dWarnick & Jeffsll2008l) 



2k h T iso B w ff R s w 



10- 



5 FM P w 



H 



R;. 



(12) 



where R; so is the correlation matrix of the array out- 
put voltages due to an isotropic external noise field at 
temperature T- lso (i.e., the correlation matrix of the ar- 
ray noise output with lossless antenna and noiseless 
receivers and with the antenna in an isotropic thermal 
environment), B is the system noise equivalent band- 
width, and R s is the array output voltage correlation 
matrix due to a signal of interest. A method for mea- 
suring R iso experimentally is described in Section l5Tl 



4.2. Noise Models 

The beam equivalen t system noise temperature is 
dWarnick & Jeffsll2008l) 
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The following sections present analytical and numeri- 
cal models used to simulate the noise covariance ma- 
trices introduced in Eq. (fj), used in Eq. ( fT3l and illus- 
trated in Fig. [5] 

4.2.1. Receiver Noise 

Because antenna elements of a PAF are closely 
spaced, mutual coupling between elements is strong. 
This influences the equivalent receiver noise tempera- 
ture through coupling of noise emitted by the input 
po rts of low-noise amplifier s (LNAs), as described 
by IWarnick & Jensen! (|2005). To account for these 
effects, the receiver noise covariance matrix R rec 
is computed in the s equel according to the model 
dWarnick et alj|2009bl) 



R rc 

where Z 
fined as 
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is the array impedance matrix, and Q is de 



T — T- 

± rcc — isc 



w H R„ 



w ff (R iso + R oss )w 



= T n 



(15) 

where Z roc is the load impedance matrix. V n , I n , 
and Y c are diagonal matrices of LNA noise parame- 
ters (noise voltage densities, noise current densities, 
and correl ation admittances, res pectively). As de- 
scribed by IWarnick et al. I d2009bl) . these terms depend 
on the minimum achievable (at perfect impedance 
match) equivalent temperature T m ; n associated with 
each LNA. In the absence of mutual coupling, R roc is 
diagonal, and each LNA can be noise matched to the 
array ports such that for all beams 



(16) 



where we have neglected downstream noise in the re- 
ceiver chains. When LNA noise couples back through 
array elements to neighboring closely packed anten- 
nas, off diagonal terms in R rcc are non-zero, and T rec 
is increased by a mutual coupling noise penalty as seen 
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in Section 15.11 Because the equivalent receiver noise 
temperature is referred to the reference plan at the ar- 
ray element ports, it includes the effects of mismatch 
loss at the interface between the array and LNAs. 

4.2.2. Spillover Noise 

This section presents a numerical model for spillover 
noise. The PAF beamformed illumination response 
pattern extends beyond the edge of the reflective sur- 
face, collecting undesired signals arriving from the 
spillover region. With a ground plane backing the ar- 
ray, backlobes are relatively small, so the spillover 
region is assumed to extend from the dish edge to the 
plane of the array. 

For a given reflector tipping angle, the spillover re- 
gion includes both sky noise and relatively high tem- 
perature (e.g., 280 K) ground. To account for this tem- 
perature distribution, we have modeled spillover noise 
as a dense grid of independent point sources with ap- 
proximately uniform angular spacing on a spherical 
annulus from the perspective of the PAF. The center 
is obscured by the primary reflector. As the dish tips, a 
portion of the annulus rises above the extended horizon 
plane, corresponding to the part of the spillover illu- 
mination pattern observing cold sky rather than warm 
ground. 

The spillover noise covariance matrix is approxi- 
mately 



R 



l6k b B 1 



J2T l Qa l a?Q H a l (17) 



where Iq is the element excitation input current, r\ is 
the intrinsic impedance of space, Ti is the noise tem- 
perature associated with the ith grid point, a; is the 
array response vector due to a unit amplitude source 
at the ith spillover noise grid position, and Ui is the 
solid angle of the corresponding sky patch. For grid 
points below the horizon, Ti = 280 K, representing 
warm ground. Above the horizon correspond, the sky 
noise model Ti = T atm (0i) + T cmh + T gb is used, 
where T atm (9) is the elevation dependent atmospheric 
noise model developed below in Eq. ( TT9b and 9i is the 
zenith angle. 

4.2.3. Sky Noise 

Atmospheric noise, CMB, and galactic background 
(GB) seen though the beamformer main lobe in the ob- 
servation pointing direction all contribute to sky noise. 



Outside the main beam, these noise sources are attenu- 
ated by the telescope's sidelobe pattern and can be ne- 
glected in the total system noise model. Atmospheric 
noise increases as the dish is tipped toward the hori- 
zon, while CMB and GB noise are modeled as a con- 
stant at all elevations. As will be shown in Section loTTl 
sky noise becomes the dominant source as the antenna 
boresight approaches the horizon. 

An extended thermal noise source larger than the 
main beam can be extended to an isotropic (i.e., full 
sphere) noise distribution with only a small perturba- 
tion to the antenna response. We can therefore use the 
approximate model to represent sky noise: 
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The isotropic noise correlation matrix R; so can be 
computed from the array element pattern overla p 
integral matrix defined by IWarnick & Jeffs d2006l) . 
Though the SOI is also seen in the main lobe, its con- 
tribution is contained in R s while R s k y includes only 
noise terms. 

For atmospheric noiase, we use a modified plane- 
parallel atmosphere model. At a given elevation angle, 
latm is proportional to the line-of-sight thickness of 
the atmosphere, which for simplicity is assumed to be 
a solid slab of uniform thickness on a flat earth surface. 
The path length through the atmosphere increases with 
depression angle 9 according to d(Q) = d sec((9), 
where do is the distance corresponding to the zenith 
direction. Since T a tm seen in the beam main lobe is 
approximately proportional to the corresponding d(9), 
we have 



{ T o . atm sec{0) 0<9<80° 

Tatm(#) = < To,atmSec(80°) 

[ +1.3(0-80°) 80° < 9 < 90° 

(19) 

where Xo ja t m = 2K is the temperature at zenith. Note 
that for depression angles greater than 80° a co rrection 



is inc luded to avoid the infinity at the horizon dRoddv 



20061) . We assume an isotropic brightness temperature 
distribution for cosmic background and galactic back- 
ground noise, with T cm b + T g b = 3K. 

4.2.4. Loss Noise 

Resistive losses in the antenna, cables, and connec- 
tors ahead of the LNA introduce a noise source that is 
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difficult to model or measure accurately. We adopt the 
approximate diagonal covariance matrix model 
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In the 19 element PAF, the dominant source of loss 
noise is a short length of coaxial line feeding each ar- 
ray element. From measurements of the cable loss, we 
have obtained the estimate Ti oss ~ 5K. 

5. Experimental Results 

5.1. Noise and Sensitivity Measurements 

The isotropic noise correlation matrix used to com- 
pute beam aperture efficiencies and system noise tem- 
peratures was measured using the ground shield facil- 
ity shown in Fig. |2] and described in Section The 
array output correlation matrix R co id for an isotropic 
cold source (sky) and Rhot with the array aperture 
covered by microwave absorber at ambient tempera- 
ture were acquired for the PAF system at 1600 MHz. 
Using an array generalization of the y-factor tech- 
nique, the isotropic noise response of the array can be 
obtained from Eq. ( TT31 in the form 



Rif 



?hot + T co \ti 



(R 



hot 



Rcold) (22) 



With the PAF mounted on the reflector, the ar- 
ray response correlation matrices for on-source, R s .i, 
and off-source pointings, R n , were obtained for each 
beam steering direction fij using calibration data as 
described in Section 13.21 These are combined with 
Rao and the known calibrator flux density F s to com- 
pute estimates of beam sensitivities, aperture efficien- 
cies and system temperature using ( [Tol l, ( TTZb . and JT31 
respectively. 

To validate the experimental measurements, a finite 
element method (FEM) numerical model of the array 
was created using HFSS (Ansoft Corp.). The ground 
plane and the dipole elements are modeled as perfect 
electric conductors, and a loss term is added to account 
for the measured loss of the coaxial cable from the 
dipoles to the element terminal connectors. 

Table Q] shows measured and modeled results for 
the 19 element prototype array for the sensitivity, aper- 
ture efficiency, and system temperature for the center 



Table 1 : Measured and modeled peak beam sensitivity, 
system temperature, and aperture efficiency for the 19 
element prototype dipole array. 

Center Element Formed Beam Model 



Sensitivity 2 m 2 /K 
T Bys 101 K 

?7ap 64% 



Table 2: System noise budget. 



3.3 irr7K 
66 K 

69% 



3.7m 2 /K 
69 K 
81% 





Measured 


Modeled 


LNA T min 


33 K 


33 K 


Mutual coupling 


20 K 


23 K 


Spillover 


5K 


5K 


Tsky 


3K 


3K 




5K 


5K 


Tsys 


66 K 


69 K 



element beam (i.e., w T = [10 ... 0]) and the full- 
array formed beam with peak sensitivity. The 12% 
discrepancy between measured and modeled aperture 
efficiency is expected, given that blockage and feed 
support scattering are not modeled and that numerical 
radiation patterns for any type of fe ed typically ov er- 
estimate efficiency by roug hly 10% dMurphvll989l) . 

The system noise temperature for a formed beam 
was estimated using Eq. ( fT3l expressed in the form 



sys 



^?rad^i£ 



w g R off w 

W^RisoW 



(23) 



The isotropic noise response essentially provides the 
available receiver gain w H Ri so w/(7y ra< ifcbTi so i?) for 
the formed beam, allowing the output noise power to 
be expressed as an equivalent sky temperature. This 
array characterization technique is particulary conve- 
nient for phased arrays, but is influenced by nonunifor- 
mities in the hot and cold noise sources and receiver 
gain variations when the PAF system is moved from 
the warm absorber/cold sky facility to the reflector. 
With a horizon model to determine the sky brightness 
temperature nonuniformity and comparing measured 
results to numerical models, the temperatures obtained 
using this technique appeared to be accurate to within 
roughly 10%. 

The system noise budget for a formed beam is 
shown in Table|2]for the 19 element array. Noise con- 
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tributions due to mutual coupling and spillover were 
not separately measured, but were estimated by car- 
rying over the modeled spillover noise and adjusting 
the mutual coupling contribution to make up the total 
measured T sys . The dominant source of loss was the 
coaxial feed for each dipole. 

The mutual coupling term in Table [2] is a contribu- 
tion to the equivalent receiver noise temperature T rec 
caused by impedance mismatches between the LNAs 
and the array. The total receiver noise temperature for 
the formed beam is T rcc = 53 K. The optimum source 
impedances of the front end amplifiers were matched 
to the isolated input impedances of the array elements, 
but due to mutual coupling between array elements the 
effective impedance presented to each front end am- 
plifier is an active impedance that is different from 
the e lement self impedance (IE. E. M. Woestenburg 
The resulting mismatches between the opti- 



2005) 



mum source impedance parameter of the LNAs and the 
active impedances presented by the array to the ampli- 
fiers led to an increase in the equivalent receiver noise. 
The active impedances are beamformer-dependent, 
which means that the equivalent receiver noise for a 
phased array varies with respect to the beam scanning 
angle. 

Improving the noise match between the front end 
amplifiers and the array active impedances to reduce 
the mutual coupling contribution to the system noise 
is a major focus of current work in PAF development. 
For communications appli cations, decoupling net- 
works have been explored (IWarnick & Jensen! 120071) 
but are likely too lossy and narrowband for astronom- 
ical instruments. More promising approaches include 
a noise matching condition that is opti mal in an aver- 



age se nse over the array field of view (IWarnick et al 
l2009albh . and the design of array antennas with ac- 
tive impedances that remain as close as possible to the 
LNA optimum source impedance (e.g., 50 il) over the 
array field of view and operating bandwidth. Results 
on the expected performance improvement with active 
impedance matching are given in Sec. 16. 21 

The values given above are for a single beam 
steered near the reflector boresight direction. All fig- 
ures of merit, including aperture efficiency and sys- 
tem temperature as well as sensitivity, vary across the 
PAF field of view. The antenna figures of merit can 
be mapped using data acquired for a grid of calibra- 
tor source pointings. For each beam, an off-source 
pointing and an on-source pointing with the reflector 
steered so that the beam is centered on the source were 
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Fig. 6. — (a) Measured beam sensitivity map (m 2 /K) 
for the 19 element dipole array on the 20-Meter Tele- 
scope (f/D = 0.43). Each pixel in the image cor- 
responds to the measured sensitivity of one formed 
beam, (b) Measured and modeled beam sensitivity for 
an elevation cut through the PAF field of view. Mea- 
sured sensitivities were obtained using a bright calibra- 
tor source on a grid of reflector steering directions with 
0.1 degree spacing. The model employed HFSS for 
the array, physical optics for reflector scattering, and a 
microwave network model for the receiver chains and 
beamformer. The half-power beamwidth (HPBW) is 
0.7 degrees. 



used to compute SNR, which can be converted to sen- 
sitivity using the known source flux density. A map 
of beam sensitivities is shown in Fig.[6{a). In this and 
all images in Section IBT21 the horizontal axis is given 
in "cross elevation" rather than azimuth, in order to 
reduce coordinate system projection distortion on the 
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sky grid. At the image midpoint, cross elevation de- 
scribes an arc perpendicular to the arc running from 
horizon to zenith. This provides an undistorted, uni- 
form pixel size rectilinear grid projection on the sky. 
Near the horizon, cross elevation is equivalent to az- 
imuth. For an image centered at zenith, both elevation 
and cross elevation axes follow arcs from zenith to 
horizon. 

Figure |6jb) shows a slice through the sensitivity 
map with a comparison to the model results. The shape 
of these sensitivity maps depends on many factors, in- 
cluding the size of the PAF, the number of elements 
used to form a beam (in this case, all 19 elements 
were used for each beam), and the location of the fo- 
cal spot in relation to the array elements. Based on a 1 
dB tolerance for sensitivity loss for steered beams and 
a raster spacing with beam overlap at the half-power 
point, the measured field of view is 2.5 HPBW in di- 
ameter. These results show good performance for the 
prototype PAF and demonstrate that scientifically use- 
ful sensitivities and efficiencies can be obtained. 

5.2. Radio Camera Imaging 

The primary motivation for PAF telescope instru- 
ments is achieving wide fields of view with multi- 
ple, electronically steered beams forming a radio im- 
age with a single dish pointing. As part of the July 
2008 experimental campaign on the Green Bank 20- 
Meter Telescope we collected a number of calibration 
data grids, used these to compute simultaneous beam- 
former weights using the maximum sensitivity beam- 
former method of ([8]), and observed several astronom- 
ical radio sources with the PAF operating as a radio 
camera. Calibration grids were either 33 x 33 or 65 
x 65 pixels in size, with inter-pixel separation of 0.1°, 
where each pixel represents a calibration pointing di- 
rection as defined in Section 13.21 For diagnostic pur- 
poses, the calibration grids were larger than the PAF 
field of view, so a subset of the grids was used to form 
images. One beam was formed per calibration pixel. 
When observing a region larger than the PAF field of 
view, the reflector was physically steered over a grid 
of pointings with approximately 1.0° spacing in order 
to form a mosaic of PAF images. 

Receiving patterns for several of the resulting 
beams are shown in Fig. [7] The peak sidelobe level 
was better than 10 dB for beams steered up to 0.6°, 
which is near the edge of the PAF field of view as 
defined by a 1 dB sensitivity loss. 
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Fig. 7. — Measured beam receiving patterns (dB rela- 
tive to peak). The 1.2° beam is beyond the PAF field 
of view as defined by a 1 dB sensitivity loss. 

Figure[8]presents examples of a single pointing im- 
age and an image mosaic. The observation is over a 
450 kHz band centered on the 1612 MHz OH line. The 
continuum source 3C295 seen in Fig. |8(a)| has a flux 
density of 21 Jy at 1400 MHz. The integration time 
was 60 seconds. The radio camera grid is 20 x 20 pix- 
els, consisting of a total of 400 simultaneously formed 
and electronically steered beams. The 3x3 mosaic 
of Fig. |8(b)| shows OH source W49N and the nearby 
3C397. Both images were oversampled with more 
beams than necessary, but they illustrate the fine-scale 
radio imaging possible. The number of pixels (beams) 
that can be formed from a data set is limited only by 
the time required to collect calibrations and the avail- 
able computational capacity, since additional beams in 
post-correlation processing do not require more array 
elements or data samples. A practical radio camera 
would likely form beams separated by on the order of 
half of the HPBW crossover distance. 

Fig.|9]presents a mosaic image of a more complex 
source distribution, the Cygnus X region, observed at 
1600 MHz with a 5 x 5 mosaic for 25 total reflec- 
tor pointings. The inset circle indicates the HPBW 
beamwidth for a single pixel beam. As a compari- 
son, the right image is from the Canadian Galactic 
Plane Survey image, but blurred by convolution with 
the equivalent beam pattern of the 20 m telescope to 
match resolution scales. We expect that the image ar- 
tifacts caused by discontinuities at mosaic tile bound- 
aries could be eliminated with more sophisticated pro- 
cessing. 
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Fig. 9. — Left: Cygnus X region at 1600 MHz. 5x5 mosaic of images using the 19-element prototype PAF on the 
Green Bank 2 0-Meter Telescope. The circle indicates the half-power beamwidth. Right: Canadian Galactic Plane 
Survey image (I A. R. Taylor et al.l 120031) convolved to the 20-Meter effective beamwidth. The center of the map is 
approximately 20 h 44 m , +42° (J2000) with north to the upper left. 



The Cygnus X radio camera image contains ap- 
proximately 3000 pixels. A more practical coarse grid 
spacing of about half the HPBW would require about 
600 pixels. A single horn feed would require 600 
pointings (one for each pixel) to form such an image, 
compared to 25 (one for each mosaic tile) for the ra- 
dio camera. Thus for equal integration times per pixel, 
this radio camera provides an imaging speed up of 24 
times. 

6. Elevation Dependent Sensitivity Optimization 

PAFs introduce the possibility of adapting formed 
beam patterns to maximize sensitivity as the noise en- 
vironment changes. Since periodic noise measure- 
ments are required to form images, the information re- 
quired to optimize beamformer weights periodically is 
already available as a byproduct of the observation. As 
the reflector pointing angle moves away from zenith, 
from the perspective of the feed the sky and spillover 
brightness temperature distributions change. In this 
section, we study the performance benefits of adap- 
tation of the beam patterns to optimize sensitivity to 
account for changes in the sky and spillover noise dis- 
tributions. 



6.1. Simulation Results 

In order to analyze the separate effects of con- 
stituent noise terms in T ays , the approach is to simulate 
expressions for R rec , R sp and R s k y at each elevation 
angle, and then compute the corresponding noise tem- 
peratures as the beamformer weights change. In or- 
der to provide an upper bound on achievable PAF per- 
formance metrics, where simulations are compared to 
a horn feed, the model incorporates active impedance 
matching (see Sec. 15. Q at the LNA-antenna interface 
to minimize effects of coupling between array ele- 
ments. This lowers the mutual coupling noise com- 
ponent of T lcc from the 23 K modeled value in Ta- 
ble|2]to 3-4 K (since active impedances depend on the 
beam steering direction, the LNAs cannot be perfectly 
matched to the array for all beam steering directions, 
and the mutual coupling component of the receiver 
noise varies over the PAF field of view). The diffi- 
culty with achieving these improvements is that active 
impedances depends on w and therefore requires a dif- 
ferent matching network for each steered beam. We 
are cu rrently developing an a rray based on the meth- 
ods of Warnick et al. ( 2009bl) which provides the best 
average match for all beams within the radio camera 
field of view and for all frequencies within the use- 
ful system bandwidth. Preliminary results suggest that 
this approach realizes most of the available improve- 
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(a) Single reflector pointing, multibeam image of 3C295. 
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(b) Mosaic image of the W49N region. 



Fig. 8. — Radio camera image examples, (a): Single 
reflector pointing, 20 x 20 pixels, 400 total simulta- 
neous beams, (b): 3x3 mosaic of PAF images ob- 
tained with nine reflector pointings. Each of the nine 
mosiac tiles has 11x11 pixels from 121 simultaneous 
beams. Lower right is 3C397, and upper left is 3C398 
(W49N). All images are in units of Jy. 



ment in T sys and sensitivity. The amplifier minimum 
noise temperature parameter T m j n was 33 K and a con- 
stant T cm b + T g b = 3 K was assumed at all elevations. 

The minimum variance distortionless response 
(MVDR) beamformer described in Section 13.31 was 
used to compute a new at each elevation to adapt 
to changes in R n . For this scenario, MVDR is es- 
sentially identical to the maximum sensitivity beam- 
former. For comparison, results are also given for a 
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(a) Mutually coupled LNA (T rcc ) and spillover (T sp ) noise 
temperatures. 
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(b) System noise temperature, T aya . 



Fig. 10. — Noise temperature variations as a function 
of dish elevation. Active impedance matching is as- 
sumed for the array. 



fixed beamformer with elevation-independent weights 
obtained from the MVDR solution computed for the 
zenith pointing. 

For the fixed beamformer curves in Fig. [TO] there 
is no fluctuation in T rcc because both the beamformer 
weights and the input receiver noise power are inde- 
pendent of the pointing direction. A decrease in T sp 
occurs as more of the pattern sidelobes in the spillover 
region are directed away from warm ground toward a 
cooler sky. T a t m (not shown) varies according to the 
secant curve described in Section 14.2.31 increasing as 
the dish elevation angle approaches the horizon and is 
identical for both the fully adaptive and fixed-adaptive 
beamformers. At Zenith T atm = 2K while at an ele- 
vation of 10° it rises to about 1 IK. 

For the adaptive beamformer, results are less intu- 
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(b) Sensitivity. 



Fig. 11. — The 19 element adaptive beamformer pro- 
vides an improvement in aperture efficiency and sen- 
sitivity as it responds to changes in the noise field. 
Higher sensitivity is achieved at mid elevations than 
for either the fixed PAF beamformer or horn feed. A 
37 element PAF would further increase sensitivity. Ac- 
tive impedance matching is assumed for each array. 

itive, since the system noise increases as the tipping 
angle moves away from zenith, as seen in Fig. |10(b)| 
As the reflector boresight moves closer to the horizon, 
the dominant change in the noise environment is the in- 
crease in the sky noise temperature in the antenna main 
beam. Since the beamformer optimizes sensitivity, the 
increase in sky noise temperature leads to an increase 
in the beam aperture efficiency. At zenith, the optimal 
beam corresponds to a reflector illumination pattern 
with a higher taper at the reflector edge, lower aper- 
ture efficiency, and higher spillover efficiency. Near 
the horizon, increased sky noise means that the opti- 
mal beam corresponds to an illumination pattern with 



a lower taper and higher aperture efficiency. This ef- 
fect is apparent in Fig. |1 l(a)| 

Figures [T0(blHrT(b)| compare PAF results for a sin- 
gle horn feed system. The horn is modeled using the 
standard cos q (<fi) illumination pattern model, where <f> 
is the angle, from the perspective of the feed, along the 
dish surface relative to the boresight axis. The value of 
q controls illumination taper and aperture efficiency, 
and is chosen to maximize sensitivity at zenith. This is 
a realistic approximate model for a typical single-pixel 
feed. For a 20 m reflector with f/D = 0.43, optimal 
sensitivity is obtained with q = 4.6, which produces 
a — 14 dB illumination taper at the dish edge. This 
comparison shows that the 19 element array, steered 
to boresight, with adaptive beamforming is capable of 
achieving a sensitivity at least as high as the typical 
single horn feed system. As noted above, to make a 
fair comparison to a single feed, for this simulation the 
LNAs are active impedance matched to the array to re- 
duce the mutual coupling noise penalty. 

Figure [Tl(b)| also includes a simulation for a 37 el- 
ement array to show sensitivity improvements that are 
possible with a larger feed. The 37 element PAF uses 
the same element design and inter-element spacing as 
the 19 element array. The significantly higher sensitiv- 
ity compared with both the standard horn feed and the 
19 element array is due to its larger aperture and the 
ability to finely control current phase and amplitude 
distribution across the array to achieve a more nearly 
optimal illumination pattern. This would be difficult to 
accomplish with standard horn feed technology, even 
with an equivalent aperture size. The larger array also 
permits beams to be steered off- axis with a smaller 
sensitivity penalty ( Waldronll2008 ). 



6.2. Experimental Results 

The simulated dish tipping scenario of the previous 
section was applied to a real data acquisition experi- 
ment with the 19 element PAF test platform on the 20m 
dish at the NRAO observatory in Green Bank, West 
Virginia. The following results use data with no bright 
sources in the array field of view collected during the 
July 2008 experiment. The observation band was cen- 
tered at 1612 MHz with a bandwidth of 450 kHz. Ar- 
ray data was recorded for telescope pointings ranging 
from zenith to the horizon in 2.5-degree increments. 

MVDR beamformer weight updates were imple- 
mented as described in Section l3~3l using to com- 
pute a new w* 3 ') at each elevation for the adaptive 
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Fig. 12. — Real data results for 19 element adap- 
tive beamforming while tipping the Green Bank 20- 
Telescope from zenith to the horizon. Simulated sen- 
sitivity is higher primarily due to scattering and block- 
age effects neglected in the model. The close curve 
shape match between simulated and measured results 
validates the simulation model. The adaptive beam- 
former exhibits improved sensitivity in the mid ele- 
vations. The physical array and simulation both used 
element-wise self impedance matching. 

beamformer. Calculations for sensitivity, aperture ef- 
ficiency, and system temperature were formed using 
elevation independent estimates of R S i and R n from 
calibration data for the boresight beam, and the eleva- 
tion dependent w (j) in ([TOj, (O and fil3l . 

Experimental results for T sys and sensitivity are 
shown in Fig.[T2l along with simulations for compar- 
ison. Simulated curves are based on the full noise 
model described in Section 14.21 Since our active 
matched array is still under development, the 19 el- 



ement test platform array and corresponding sim- 
ulations reported in this section used element-wise 
self impedance matching (i.e., the LNAs are noise 
matched to the element self impedances, not the active 
impedances). Comparison of Figs. |12(a)| and |10(b)| 
illustrates that active impedance matching (see Sec. 
15.11 ) has the potential for this array of reducing T sys 
by 23 K. Ohmic losses in Fig. 10(b)| are assumed to 
be negligible, whereas Fig. |12(a)| includes 4K noise 
due to loss. Figures |12(b)| and 1 l(b)| show that ac- 
tive impedance matching can increase sensitivity by a 
factor of 1.67 (2.2 dB). 

There is significant agreement in Fig. [T2l between 
experimental real data and simulated results, which 
serves as a validation for the proposed noise models. 
In particular, the increase in T sys , ?7 ap , and the modest 
improvement in sensitivity at lower elevations for the 
adaptive beamformer with respect to the fixed beam- 
former predicted by the model simulations has been 
verified with real data. Computed real data aperture ef- 
ficiency (not shown) also exhibited an excellent match 
with the model simulations. The good match in this ex- 
ample application suggests that the model may be used 
with confidence in a variety of other signal scenarios. 

In some observing scenarios it may not be practi- 
cal to obtain an SOI-free RW estimate to update w^^ 
continuously with elevation changes. In our example 
we took pains to find a relatively source-free patch of 
sky for the elevation scan. Since elevation dependence 
of sky and spillover noise is based primarily on the 
geometry and not on transient sources, however, it is 
possible to pre-compute a set of optimal, elevation in- 
dexed beamformers from calibration data. These could 
then be called up at a later time (over several days) 
from a look-up table to achieve the improved mid- 
elevation sensitivity shown above. 

The degradation rate of calibration data due to 
phase and amplitude drift in the receiver electronics 
and other effects is currently under study, as are tech- 
niques to refresh calibration data using a small num- 
ber of bright source observations or external electronic 
calibration sources. 

7. RFI Mitigation 



As contemporary science goals increasingly require 
observing sources outside the traditional protected 
spectrum bands, a critical need is developing to deal 
with ubiquitous man-made inter fering signals such as 
satellite downlink transmissions ( S.W. Ellingsonl2001 
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air n avigation aids dZhang et al 
20051) . wireless communications (Leshem et al 



1994), radar 



Poulsen et alj 20051 Combrinck et alj 
systems foong et al.l 20051: jjeffs et al.1 l2006t iFisher 
200ll: lEllingson & Hampsonl2003): | Zhang et all2003l) 



20051: 



Fisher et al.. 
19991: 



Boo nstra et al. l2000h . and digital television broad- 



casts. Even locating instruments in undeveloped ar- 
eas with regulatory protection does not avoid many 
man-made sources such as satellite downlinks. 

The PAF, as illustrated in Fig. |4j offers a promising 
new approach that exploits the spatial structure of the 
interfering signa l to track and remo ve it without having 
to discard data ( Jeffs et al. 2008bi) . With spatial can- 
celing, even interferers that entirely overlap the signal 
of interest (SOI) spectrum may be mitigated. It is an- 
ticipated that using this technique, data collection may 
now be possible where previously interference was too 
dominant to permit serious scientific observations. 

To demonstrate feasibility of adaptive RFI miti- 
gation with the 19 element array on the 20m dish, 
a local, moving, man-made interference source was 
introduced. This involved an RF function generator 
and hand-held antenna, while walking a circuit located 
about one kilometer from the telescope. The interfer- 
ing signal was FM-modulated with a 200 kHz band- 
width overlapping and masking the W30H spectral 
line seen at the array at 1665 MHz. RFI was removed 



using t he subspace projec tion algorithm (iLeshem et al 



2000; Ueffs et al.l |2008bJ). This method adapts the 



beamformer weights to the changing interference envi- 
ronment, so as to place a spatial null in the interference 
direction. Beamformer weights are updated every STI 
according to 



P x = 



P J 
I 



U in t.U 



R w [U int |U n ] = [U int |U n ]A 



(24) 



where j is the STI index, sample covariance Rw is 
computed as in and Wj is the nominal determinis- 
tic (non adaptive canceling) beamformer weight vector 
for the beam steered in the ith direction. Equation d24l i 
is an eigenvector decomposition of R") with eigen- 
values in diagonal A ordered in descending magnitude 
such that the dominant eigenvectors in the partition 
Ui n t span the interference subspace. The full eigen- 
vector matrix U = [Uj nt | U n ] is unitary. 

Images of the source with and without RFI miti- 
gation are shown in Fig. [13] Some small distortion 
due to residual RFI is apparent in Fig. \l3[c) after sub- 




(a) No RFI. 




(b) With RFI. 




(c) RFI canceled. 



Fig. 13. — W30H image with and without RFI. The 
color scale is equivalent antenna temperature (K). 



space projection adaptive cancelation processing, but 
the source which was completely obscured by interfer- 
ence is now clearly visible. 

8. Conclusions 

We have reported experimental results for a proto- 
type L-band PAF and have demonstrated radio camera 
imaging using in situ calibrated maximum sensitivity 
beamforming. Experimental results agree with simu- 
lations for a complete electrical and noise model for 
the PAF and reflector system. Measured beam sen- 
sitivities, aperture efficiencies and system noise tem- 
peratures achieve expectations for the 19 element, un- 
cooled array. The potential of PAFs for adaptive RFI 
mitigation was explored. These results indicate that 
there is a clear path forward to a high sensitivity, wide 
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field of view PAF with scientifically useful perfor- 
mance. 

Ongoing work includes experimental demonstra- 
tion of active impedance matching to reduce the noise 
contribution caused by array mutual coupling, studies 
of beamformer design procedures with controlled pat- 
tern shape, development of an L-band PAF and asso- 
ciated back-end signal processing for the Green Bank 
Telescope, and cryogenic arrays that achieve system 
temperatures of approximately 20 K. 

This work was funded by National Science Foun- 
dation under grant number AST 0352705. We express 
appreciation to Jay Lockman for providing examples 
of the science that will be enabled by phased array 
feeds. 

Appendix 

The following notations are used in the paper: 

1. E{A} : Expected value of random A. 

2. E m (Tl): Far-field electric field pattern at spheri- 
cal angle SI due to array element m. 

3. L a J : Floor operation, rounding toward zero. 

4. a* : Complex conjugate of a. 

5. A T , A H : Transpose and conjugate transpose 
of A. 

6. A : Sample estimate of some parameter, A. 

7. A : Calibration data set estimate of A. 

8. I : The identity matrix. 
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